EXPONENTIAL PREFIXED POLYNOMIAL EQUATIONS 



ARAN NAYEBI 

Abstract. A prefixed polynomial equation is an equation of the form P(ti, . . . ,t n ) = 0, where P is a 
polynomial whose variables t\ , . . . , t n range over the natural numbers, preceded by quantifiers over some, or 
all, of its variables. Here, we consider exponential prefixed polynomial equations (EPPEs), where variables 
can also occur as exponents. We obtain a relatively concise EPPE equivalent to the combinatorial principle 
of the Paris-Harrington theorem for pairs (which is independent of primitive recursive arithmetic), as well as 
an EPPE equivalent to Goodstein's theorem (which is independent of Peano arithmetic). Some new devices 
are used in addition to known methods for the elimination of bounded universal quantifiers for Diophantinc 
predicates. 
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1. Introduction 

A prefixed polynomial equation is an equation of the form P(t\, . . . ,t n ) = 0, where P is a polynomial 
with variables t\,...,t n that range over the natural numbers, preceded by quantifiers over some or all of 
its variables. Bovykin and De Smet [1] study the collection of all such possible expressions (terming this 
"the Atlas"), with the equivalence of relation of being "EFA-provably equivalent" on its members. Thus, 
members of the same class are prefixed polynomial expressions that are provably equivalent to one another. 
It is not difficult to obtain a prefixed polynomial representation, but the value of obtaining polynomial 
expressions is that they provide concrete examples of unprovable statements and explicit illustrations of 
deep logical phenomena. For example, the prefixed polynomial expression that Bovykin and De Smet obtain 
for l-Con(ZFC+Mahlo) implies, over /Si, all two quantifier arithmetical theorems that can be proved in 
ZFC + Mahlo cardinals. To avoid too much repetition, we refer the reader to the discussion in [1] for a 
detailed exposition as to the importance of such an Atlas. 

One example of prefixed polynomial expressions of well-known logical phenomena that Bovykin and De 
Smet obtain are representations (involving alternations of universal and existential quantifiers) of the Paris 
Harrington theorem, and the special cases of the Paris-Harrington theorem for pairs and triples. The Paris- 
Harrington theorem [15] , which states that a simple extension to the finite Ramsey theorem is not provable 
in first-order Peano arithmetic, is a seemingly natural mathematical example of incompleteness, namely: 



Date: September 12, 2012. 

2010 Mathematics Subject Classification. Primary 03B70; Secondary 05C15, 05A18. 
Key words and phrases. Diophantine equation, Finite Ramsey Theorem. 



1 



2 



ARAN NAYEBI 



Theorem 1 (PH). For all numbers e, r, and k, there exists a number M, such that for every coloring f of 
e-subsets of [M + 1] = {0, 1, . . . , M} into r colors, there is an f -homogeneous Y C [M + 1] of size at least 
min(y) + k — 1. 

Our focus will be on the Paris-Harrington theorems for pairs, whose combinatorial principle will be referred 
to here on out as PH 2 : 

Theorem 2 (PH 2 ). For every number k, there exists a number M such that for every coloring f of 2- sub sets 
of [M + 1] into r colors, there is a f -homogeneous Y C [M + 1] of size at least min(Y) + k — 1. 

For PH 2 and r > 2, (where r is the only free variable which represents the colors), the Bovykin-De Smet 
prefixed polynomial representation |T[ Theorem 2] is: 

Vfc 3M Vab BcdAX Vxy 3BCF Vfg Behilnpq 

[x-(y + B-x)-(A + k + B-y). ((((/ - A) 2 (g - if) ■ ((/ - B) 2 + (g - xf) 
' ' ' ' ' ■ ((/ - Cf + (g - yf) -h-1)- {{dgi + i-c+f) 2 + (/ + h- dgf) + (B + l + l-C) 2 

+ (C + n- Mf + (F + e-b(B + C 2 )) 2 + (bp(B + C 2 ) + p - a + F) 2 + {(F - X) 2 - qr) 2 )]. 

Although in the case of PH 2 , the prefixed polynomial expression covers only a few lines, the challenge 
comes when transforming this polynomial from its n[! form to its EFA-provably equivalent n° form. Here, 
bounding the universal quantifiers and then eliminating them introduces a drastic increase in the number 
of variables of the original prefixed polynomial representation, to the point that the resulting n° form is 
too long to be practical to write. The transformation of formulas containing bounded universal quantifiers 
into equivalent formulas containing only existential quantifiers in the theory of Diophantine equations is 
a powerful technique which has many applications, such as showing in a straightforward manner that the 
set of primes is Diophantine, constructing a universal Diophantine equation, or demonstrating that many 
famous problems can be reformulated in terms of the unsolvability of a particular Diophantine equation 
(since many of these problems such as Goldbach's conjecture, the Riemann hypothesis, and the four color 
theorem can be formulated in the form Vn P(n), where P is a decidable property over natural numbers). 
However, naive attempts to obtain a Diophantine representation (namely, a direct application of the results 
of Davis, Putnam, and Robinson [2] and Matiyasevich [8], and possibly with some slight modifications but 
with no drastic tricks) for PH 2 yields unwriteable representations. 

We discuss the methods that are used to eliminate the bounded universal quantifier and present several 
ways of conserving the large number of variables typically introduced by this process in order to obtain the 
following result: 

Theorem 3 (Unprovability by primitive recursive means). There exists an exponential Diophantine equation 
Ei(k, M, a, b, r, ti, . . . , ii3s) o,nd a Diophantine equation D\(k, M, a, 6, r,ti, . . . , £347) (both with k, M, a, 6, and 
r as parameters) such that for every r > 2 

Ei(k, M, a, 6, r, ii, . . . , t 138 ) = 

has a solution in natural numbers t\, . . . , t\$% and 

Dx(k, M,a,b,r,ti,..., t 347 ) = 

has a solution in natural numbers ti,...,t347 is equivalent to the combinatorial principle of the Paris- 
Harrington theorem for pairs, equivalent to the 1-consistency of /Si and thus not provable in /Ei (but 
provable in IT>2). 

We will prove Theorem |3] in §3.1 and provide an explicit representation of Ex(k, M, a, b, r, t±, ti^s). In 
§3.2, we consider unprovability in Peano arithmetic and obtain an explicit EPPE equivalent to Goodstein's 
theorem, obtaining the following result: 

Theorem 4 (Unprovability in Peano arithmetic). There is a 181 variable exponential prefixed polynomial 
equation equivalent to Goodstein's theorem, unprovable in Peano arithmetic. 
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Our results offer insight in preserving a writeable representation when we reduce the quantifier complexity 
of the original prefixed polynomial representations of Bovykin and De Smet [1 , therefore not restricting us 
to using alternations of universal and existential quantifiers in order to explicitly illustrate deep logical 
phenomena. Not only that, but our consideration of exponential prefixed polynomial equations allows one 
to obtain short representations of "natural" independent statements for which exponentiation is inherent in 
the formulation, such as Goodstein's theorem. 

2. Elimination of the bounded universal quantifier 

There are several methods of transforming formulas with bounded universal quantifiers to those having 
only existentially bound variables in the language of Diophantine predicates. The most well-known is the 
Bounded Quantifier Theorem of Davis, Putnam, and Robinson [2], which uses the Chinese remainder theorem 
to establish this equivalence. Bounded universal quantifiers can also be eliminated by way of Turing machines 
(presented in detail in Chapter 6.1 of Matiyasevich P3]), a rather immediate consequence of Matiyasevich's 
direct method in |12j of simulating Turing machines by Diophantine equations. Finally, a third method of 
going about this elimination involves summations of generalized geometric progressions, based on a technique 
first proposed by Matiyasevich in [13] and presented for the first time in Chapter 6.3 of Matiyasevich's book 
[13] . There are advantages and disadvantages to each method. The method via Turing machines, though 
constructive, is rather roundabout. With the method involving summations, although straightforward, 
it becomes impractical to extract the resultant Diophantine equation from the heavy use of generalized 
geometric progressions after the elimination of the bounded universal quantifier when the expression obtained 
prior to this elimination is of even moderate size. Thus, this is our motivation for preferring the method of 
Davis, Putnam, and Robinson via the Chinese remainder theorem because it is a straightforward number- 
theoretic method which produces, in most cases, a visualizable Diophantine equation. However, as we 
will explicitly demonstrate below, this does not mean that the representation can be, practically speaking, 
explicitly written down since the downside is the drastic number of variables introduced. 

For example, for the following: 

(2.1) \fy < b 3xi, ■ ■ ■,x m [G(a,y,Xi, . . .,x m ) = 0], 

where a represents the parametcr(s) of the polynomial G, the Chinese remainder theorem method results in 
the following system of Diophantine conditions solvable in the unknowns q, w, zq, . . . , z m provided that (|2.ip 
holdfi 



G{a,ZQ,Zi,...,z m ) = mod 

(2.2) 

b\{b + w + B{a,b,w))\ \q + l, 
q\ , { zi 



if 



(2.3) 



where the polynomial B(a, b, w) is obtained from G(a, y, xi, . . . , x m ) by changing the signs of all its negative 
coefficients and systematically replacing y by b and Xi, . . . , x m by w. 

Now, the major contributing factors to the increase in the number of variables are the representations of 
the factorial and binomial coefficient. Let T denote the number of variables in the Diophantine representation 
of the exponential function, then the Diophantine representations of the factorial and binomial coefficient 



^This is based on the original construction of the Bounded Quantifier Theorem of Davis, Putnam, and Robinson [2], with 
some minor modifications by Matiyasevich [141 Ch. 6.2], There have been more drastic modifications to the Bounded Quantifier 
Theorem, namely by Matiyasevich |9l |10| and by Hirose and Iida 5 , and we incorporate some of these results in our presentation. 
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as presented by Davis, Matiyasevich, and Robinson [31 §1] involve 10 + 5(r + 1) variables and 6 + 3(r + 1) 
variables, respectively. And the Diophantine representations of the factorial and binomial coefficient as later 
presented by Matiyasevich [14l Chapter 3.4] involve 10 + 6(r + 1) variables and 5 + 4(r + 1) variables, 
respectively. If we used an economical (with respect to the number of variables) Diophantine representation 
of the exponential function, for instance, the result obtained by Matiyasevich and Robinson with only 
five variables, then the total number of variables for the factorial and binomial coefficient presented in [3J 
would be 40 variables and 24 variables, respectively; and the results obtained in [T3] would be 46 variables 
and 29 variables, respectively. Now, note that the representation of the exponential function in five variables 
obtained in results in a polynomial of high degree and is a somewhat unruly expression (though obviously 
writeable) , for both the binomial coefficient and the factorial. Furthermore, it is important to note that (|2.3[) 
involves m binomial coefficients (I am excluding m). Thus, even if we used the least number of variables in 
representing the binomial coefficient, 25 variables (since for each binomial coefficient, we have to introduce 
a new variable yi, such that yi = ( Zl ), for each / = l,...,m), then (|2.3[) would require 25m variables to 
represent it. For m = 3, a possibly small Diophantine equation with a single universally bound variable y, 
this would mean the introduction of 75 new variables just to represent the system of conditions in (|2.3|) ! 

Thus, the first step toward mitigating the increase of variables introduced by (|2.2j) and (|2.3[) the number 
of variables involved in the Diophantine representations of the factorial and the binomial coefficient must 
be drastically reduced. Fortunately, it turns out that one can eliminate the need to show that factorial is 
Diophantine. By a result of Matiyasevich [9] [10], one can use the multiplicative version of Dirichlet's box 
principle to replace the condition in (|2 . 2|) : 

(2.4) bl(b + w + B(a,b,w))l | q+ 1 
with the sufficiently strong inequality, 

(2.5) q > b + (b + l) b+1 ((b + l) b+1 B{a, b, w)) wT " . 

As can easily be seen, (|2.4[) would require the introduction of 10(r + l) + 22 variables (using the representation 
of the factorial provided by Davis, Matiyasevich, and Robinson [3) §1]) whereas (|2.5p would only require the 
introduction of 2(r + 1) + 2 variables. Hence, even if one is content with just an exponential Diophantine 
representation, then 20 variables are already conserved, and if one would prefer a Diophantine representation 
even with an economical representation of the exponential function in only five variables, then 68 variables 
are conserved! 

However, what has not really been proposed so far is a reduction in the number of variables introduced 
by (|2.3p . the system of binomial coefficients, since that is the primary reason why so many variables are 
introduced in the elimination of the bounded universal quantifier (as demonstrated above in the case for 
m = 3). We will prove that one can represent the binomial coefficient in only 10 variables and obtain an 
explicit representation that covers less than a page at 18 variables. 

In our construction of a representation of the binomial coefficient in only 10 variables, we will rely on 
the relation-combining theorem of Matiyasevich and Robinson [TT]. The Matiyasevich- Robinson relation- 
combining theorem allows one to cheaply define certain combinations of relations than by defining each 
separately by an equation and then combining the equations. While it is economical with respect to the 
number of variables of the resultant equation, it should be noted that the relation-combining theorem is 
rather uneconomical with respect to the degre^E The theorem is stated as follows: 

Theorem 5 (Relation-combining theorem) . Let □ denote a perfect square. For all integers A\ , . . . , A q , B, C, D 

with B ^ 0, the conditions Ai — □ (i = 1, . . . , q), B | C, and D > all hold if and only if M q (A±, . . . , A q , B, C, D, n) — 

for some n, where M q is the following 2 q -fold product over all combinations of signs 

M q = Y[(B 2 n + C 2 - B 2 (2D - 1) ■ (C 2 + W q ±^/a\±^/A~ 1 W ± . . . ± v^tf/ 9 " 1 )), 



A version of the relation-combining theorem that is more efficient with respect to the degree was later worked out by 
Matiyasevich, presented as Theorem 5.1 of Jones [6]. However, this does not change the fact that the resultant polynomial can 
cover several pages, which is why we rely on more elementary techniques to have our representation of the binomial coefficient 
cover half a page, with an introduction of only eight more variables. 
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where 



We also need a result on the solutions to the Pell equation 
(2.6) x 2 - (a 2 - l)y 2 = 1. 

For a > 0, we define the pair < Xa(fi), ^a(^) > as the n-th nonnegative solution of (|2.6p . Theorem 4 of 
Matiyasevich and Robinson proves the following system of Diophantine conditions: 

Theorem 6. For A > 1, B > 0, and C > 0, C = i/ja{B) if and only if the following system of conditions is 
satisfied: 

DFI = □, F\H-C, B <C, 

D = {A 2 -1)C 2 + 1, 

E = 2(i + l)DC 2 , 

F = (A 2 - l)E 2 + 1, 

G = A + F(F — A), 

H = B + 2jC, 

I = (G 2 — l)H 2 + 1. 

From Theorems [5] and [6l we are led to the following result: 

Theorem 7. The relation y = ( n j , where n > s > 0, holds if and only if 

F(y,n,s,x,w,k,l i m,i,j,vi,v 2 ,v 3 ) = 

has a solution in the parameters y,n, and s and the variables x,w,k,l,m,i,j,Vi,V2, and V3. 
Proof. It is easy to see that 



(2.7) 



y = 

^ 3x 
<^ 3xw 



y = rem 1 



(x + iy 



y = rcm(i« + A W + 1 = 
where the partial binomial expansion 



x ) A x > An 3 

{x + iy 



Ax - An 8 > 



is defined as 



(x + iy 



71— 8 

\ - I n 



i=0 



S + I 



for n > 0, s > 0, and x > n s . Note that the reason we take w + 1 as opposed to just w in (|2.7p is to ensure 



that w + 1 > 0. By Theorem 8 of [IT], w + 1 



can be expressed as a system of Diophantine 
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conditions in three variables, so (|2.7|) becomes the following system of equations 



y = rem(w + l,x), 

C = i/> A (B), 

(M 2 - l)K 2 + 1 = □, 

(M 2 x 2 - 1)L 2 + 1 = □, 

(x - 4n s )(K 2 L 2 - 4(C - KL(w + l)) 2 ) > 0, 

M = 8n(x + w + 1) + 2, 

K = n- s + l + k(M - 1), 

L = s + l + l(Mx - 1), 

A = M(x + 1), 

B = n+1, 

C = m + B 

Note that 

y = rcm(w + 1, x) 

(y + Vi + 1 - x) 2 + (xv 2 - w - 1 + y) 2 = 0. 

Moreover, since C > B, the remaining conditions can be combined via Theorem [5] as 

M 3 (DFI, (M 2 - l)K 2 + 1, {M 2 x 2 - l)L 2 + 1,F,H — C,(x — An s ){K 2 L 2 - 4(C - KL{w + l)) 2 ), v 3 ) = 0. 

Since we can eliminate D, F, I, H, M, K, L, A, B, and C, we get that 

F(y, n, s, x, w, k, I, m, i, j, v 1 ,v 2 ,v 3 ) 
= ((y + v 1 + l- x) 2 + (xv 2 -w-l + y) 2 ) 2 

+ Ml{DFI, (M 2 - l)K 2 + 1,(MV - l)L 2 + l,F,H-C, (x - 4n s ){K 2 L 2 - 4(C - KL{w + l)) 2 ), v 3 ). 



The explicit representation of F(y, n, s, x, w, k, I, m, i,j, v±,V2, V3) is rather unruly and would cover several 
pages. Instead, we keep D,F,I,K, L,M as variables and introduce two more, W and J, where W = 
1 + (DFI) 2 + ((M 2 - l)K 2 + l) 2 + ({M 2 x 2 - l)L 2 + l) 2 and J = (x — 4n s )(K 2 L 2 - 4(C - KL(w + l)) 2 ). 
Thus, we combine the additional equations via the summing of squares technique, and the following explicit 
representation of Fi(y,n,s,x,w,k,l,m,i,j,Vi,V2,V3,D,F,I,J,K,L,M,W) is obtained, where y,n, and s 
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are parameters: 



((y + v x + 1 - x) 2 + {xv 2 -w-l + y) 2 ) 2 + {{D - ((M(x + l)) 2 - l)(m + n + l) 2 - l) 2 + (F - 4((M(x + l)) 2 - 1) 

(i + l) 2 (((Af (a + l)) 2 - l)(m + n + l) 2 + l) 2 (m + n + l) 4 - l) 2 + (/ - ((M(x + 1) + (4((M(x + l)) 2 - l)(i + l) 2 

(m + n + 1) 4 (((M(.t + l)) 2 - l)(m + n + l) 2 + l) 2 )(4((M(a; + l)) 2 - l)(i + 1) 2 (to + n + l) 4 (((M(x + l)) 2 - 1) 

(m + n+ l) 2 + l) 2 - M(a: + 1) + l)) 2 - l)(n + l + 2j(m + n+ l)) 2 - l) 2 + (A/ - 8n(a; + w + 1) - 2) 2 + (K-n + s-l 

- k(M - l)) 2 + (L - s - 1 - i(Mrs - l)) 2 + (W - 1 - (DFI) 2 - ((A/ 2 - 1)X 2 + l) 2 - ((M 2 x 2 - l)L 2 + l) 2 ) 2 

+ (J — (x — An s )(K 2 L 2 - 4(m + n + 1 - KLfra + l)) 2 )) 2 ) 2 + ((F 2 w 3 + (n + 1 + (m + n + l)(2j - l)) 2 - F 2 (2J - 1) 
((n + 1 + (m + n + l)(2j - l)) 2 + W/ 3 + (DFI) 1/2 + ((M 2 - 1)K 2 + 1) 1/2 W + ((AfV - 1)L 2 + l) 1/2 W 2 )){F 2 v 3 
+ (n + 1 + (m + n + l)(2j - l)) 2 - F 2 (2J - l)((n + 1 + (m + n + l)(2j - l)) 2 + W 3 + (DFI) 1/2 - ((A/ 2 - l)K 2 + 1) 1/2 
W + {{M 2 x 2 - l)L 2 + l) 1/2 W 2 )){F 2 v a + {n + 1 + (to + n + l)(2j - l)) 2 - F 2 (2J - l)((n + 1 + (to + n + l)(2j - l)) 2 
+ W 3 + (DFI) 1/2 - {{M 2 - l)K 2 + l) 1/2 iy - ((A/V - 1)L 2 + l) 1 / 2 W 2 ))(F 2 v 3 + (n + 1 + (to + n + l)(2j - l)) 2 

- F 2 (2J - l)((n + 1 + (TO + n + l)(2j - l)) 2 + T4^ 3 - {DFI) 1 / 2 - ((A/ 2 - l)K 2 + \f/ 2 W - {{M 2 x 2 - 1)L 2 + lf/ 2 W 2 )) 
(F 2 v 3 + (n + 1 + (to + n + l)(2j - l)) 2 - F 2 (2J - l)((n + 1 + (to + n + l)(2j - l)) 2 + W 3 - (DFI) 1/2 - ((M 2 - l)K 2 
+ 1) 1/2 W + {{M 2 x 2 - l)L 2 + l) 1/2 W 2 )){F 2 v a + (n + 1 + (to + n + l)(2j - l)) 2 - F 2 (2J - l)((n + 1 + (to + n + 1) 

(2j - l)) 2 + Ty 3 - (DFI) 1/2 + {{M 2 - 1)F 2 + l) 1/2 iy + ((A/ 2 a; 2 - l)L 2 + l) 1/2 Ty 2 ))(F 2 w 3 + (n + 1 + (m + n + 1) 
(2j - l)) 2 - F 2 (2J - l)((n + 1 + (to + n + l)(2j - l)) 2 + VK 3 - {DFI) 1/2 + ((A/ 2 - l)^ 2 + 1) 1/2 T^ 

- ((Af 2 a; 2 - 1)L 2 + l) 1 / 2 W 2 ))(F 2 v 3 + (n + 1 + (to + n + l)(2j - l)) 2 - F 2 (2J - l)((n + 1 + (to + n + l)(2j - l)) 2 
+ W 3 + (DFI) 1/2 + ((A/ 2 - l)K 2 + l) 1/2 iy - ((Af 2 a; 2 - l)L 2 + 1) 1/2 T^ 2 ))) 2 . 



Now, if one would prefer to have a writeable exponential Diophantine equation from the transformation of 
a formula with bounded universal quantifier(s) (since for certain problems obtaining a writeable exponential 
Diophantine representation is more feasible than obtaining a writeable Diophantine representation, as will be 
the case in §3), then even the representation of the binomial coefficient in 18 variables is too large, since one 
can use the usual exponential Diophantine representations of these in 5 variables and 10 variables, respectively 
(if one uses the exponential Diophantine representations presented in [T31 Chapter 3.4]). However, the 
problem one still faces is the vast number of variables introduced by (|2 . 3[) . The relation of divisibility 
requires the introduction of one new variable, and each binomial coefficient (again, I am excluding (?) from 
this) requires the introduction of six new variables. Thus, 7m new variables are needed in the exponential 
Diophantine representation of (|2.3|) . For example, in §3, we will be dealing with m — 24 and to = 31, which 
would mean for each case 168 and 217 variables are introduced, respectively. 

Therefore, a method of reducing the number of variables involved in (|2.3[) would do us well. More 
specifically, the goal is to reduce the number of variables in 



where yi = (?) and I = 1, . . . , to. 

Inspired by some tricks used in the proof of the result by Matiyasevich [13] that every Diophantine set 
has an exponential Diophantine representation with only three unknowns, the following result is obtained: 



(2.8) 




Theorem 8. 



(2.9) 




3 P q[((yi + 1)2 Z > + l) Zl = P (( yi + 1)2*T + 9 A (( yi + 1)2 2 



> Q A y x | p] . 
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Hence, since the relations ((yi + 1)2 Z ') W > q and y\ \ p introduce 2 more variables, we have a total of 4 new 
variables introduced by (|2.9[) . 

Proof. Observe that by the binomial theorem, 

(2.10) (u + l) 2 ' = pu w + q, 

where 

zi 



E 

i—w 
w — 1 

E 



If u is large enough, for instance, if 
then 



u > 2 Zl , 



u w > u w - x ¥ l 



=e : 







> q- 



Note that p and q are uniquely determined by (|2.10l) and (12. lip . 

It is easy to see that since 



E 



then 



w I 1 — ' \ i 

i—w+l 



p = [ 1 I mod u. 



So if ?/i | u, then the condition in (|2.8p is equivalent to 



2/i I V- 



Thus, 



^ 3pqu[(u + l) Zl = pu w + q A u > 2 2 ' A u w > q A yi \uAyx p]. 



The conditions u > 2 Zl , u w > q, y\ \ y, and y\ \ p introduce 4 more variables, bringing the total again, to 7 
new variables. However, we can reduce the total number of new variables introduced, namely by eliminating 
u and 2/i | u by using the equation 

(2.12) u=( yi + l)2 Zl . 

Thus, we get our result. ■ 

So, the representation of (|2.8p has been reduced from introducing 7 new variables, to introducing only 4 
new variables. Hence, already for m > 7, we see that the number of variables conserved by (12.91) supercedes 
the number of variables conserved by the strong inequality in (|2.5[) . 
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3. Exponential prefixed polynomial representations of independent statements 

3.1. Proof of Theorem[3j Bovykin and De Smet's [T] intermediate representation of Theorem[5]in prefixed 
polynomial form (involving alternations of existential and universal quantifiers) is as follows: 



Vfe 3M Va6 BcdAX \/xy 3BCF 
[(Q <xAx <y Ay <A + k-l) ^ 

(A = rcm(c. d + 1) A B — rem(c, dx + 1) A C = rem(c, dy + 1) A B < C 
AC < M + 1AF = rem(a, b(B + C 2 ) + 1) A F = X mod r)]. 



Note that eliminating the bounded quantifiers from the intermediate representation of Theorem[2]rather than 
directly from the full prefixed polynomial in (jl.ip results in smaller exponential Diophantine representation 
for PH 2 . Our notation is pretty similar to [I], but to clarify, (a, b) codes / and (c, d) codes the homogeneous 
set Y (and x and y are indices of elements coded by (c, d)). Moreover, A is the first element of Y, namely 
min(y), since Y is ordered. 

Note that we will be taking k, M, a, b, and r as parameters (thus, the Vfc and Va6 in (|3.ip pose no difficulty). 
Further, note that since a; is a natural number x > 0, but in (|3.ip . x > 0. To fix this, we simply have to 
modify B and C to B — rem(c, d(x + 1) + 1) and C = rem(c, d(y + 1) + 1). Second, since we have updated x 
and y to be x + 1 and y + 1, respectively, and x < y, then Vxy < A + k — 1 is incorrect. Instead, we bound 
x and yasVx<A + k — 3 Vy < A + k — 2. Hence, the representation we are left to deal with is: 



BcdAX Vx<A + k-3Vy<A + k-2 3BCF 
(3.2) [x < y A A = rem(c, d + 1) A B = rcm(c, d(x + 1) + 1) A C = rcm(c, d{y + 1) + 1) A B < C 

A C < M + 1 A F = rem(a, b(B + C 2 ) + 1) A F = X mod r}. 



Expanding ()3.2|) . we have 



3cdAX \/x<A + k-3Vy<A + k-2 3BCFv u ...,v 17 
[{x + vi + 1 - y) 2 + {{A + v 2 - df + ((d + l)wa - c + A) 2 ) 2 
+ ((B + v 4 - d{x + l)) 2 + ((d(x + 1) + 1)« B - c + S) 2 ) 2 
+ ((C + « 6 - d(y + l)) 2 + ((d(y + 1) + 1)«7 - c + C) 2 ) 2 
+ (B + w 8 + 1 - Cf + (C + v Q - Mf 

+ ((F + « 10 - b(B + C 2 )) 2 + ((b(B + C 2 ) + l)v n - a + F) 2 ) 2 

+ ((«12 + «14 + 1 - ^) 2 + 0U15 - F + V 12 ) 2 ) 2 

+ ((vis +v 16 + l- r) 2 + (rv 17 -X + v 13 ) 2 ) 2 + (v 12 - v 13 ) 2 = 0]. 



We can reduce the two bounded quantifiers in (|3.3[) to just one by taking advantage of the fact that if 
x < A + k — 3 and y < A + k — 2, then J(x, y) < J (A + k — 3,A + k — 2), where J is Cantor's function 
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defined for natural numbers m and n as J(m, n) — \{(m + n) 2 + 3m + n). Thus, we have 



BcdAXz Mt < z + 1 BxyBCFvi,. ..,vw 

[(22 - (2A + 2k- 5) 2 - 4A - 4k + ll) 2 + (2t - (x + y) 2 - 3a; - y) 2 

+ ((A + A; - 2 + u 18 - x) ■ (A + k - 1 + v w - y) ■ ((x + v x + 1 - y) 2 

+ ((A + v 2 - d) 2 + ((d + l)v 3 - c + A) 2 ) 2 

+ ((B + v 4 - d(x + l)) 2 + {{d{x + 1) + 1)» B - c + B) 2 ) 2 

+ ((C + w 6 - d(y + l)) 2 + ((d(y + 1) + 1)«7 - c + C) 2 ) 2 

+ (B + v 8 + 1 - C f + (C + v 9 - M) 2 

+ ((.F + Vio - KB + C 2 )) 2 + ((b(B + C 2 ) + l)« u - a + ^) 2 ) 2 
+ ((«12 +ui4 + 1 - r) 2 + (ruis - F + wi 2 ) 2 ) 2 

+ (("is + vie + 1 - O 2 + (rv 17 -X + v 13 ) 2 ) 2 + (v 12 - v 13 ) 2 )) 2 = 0]. 



The final removal of the remaining bounded universal quantifier Vt < z+1 will be done explicitly by the meth- 
ods presented in §2. Based on (|3.4|) . we then define the following polynomial P(k, M, a, 6, r, c, d, A, X, 2, z , 21, . . 
as: 



(2z - (2A + 2k- 5) 2 - 4A - 4fc + f f ) 2 + (2z - {zi + z 2 ) 2 - 3zi - z 2 ) 2 

+ ((A + fc - 2 + z 2 3 - zi) • (-4 + fc - f + z 24 - 22) • ((21 + z 6 + 1 - z 2 ) 2 

+ ((A + z 7 - d) 2 + ((d + l)z s - c + A) 2 ) 2 

+ ((z 3 + z 9 - d(z x + I)) 2 + ((d(zi + 1 ) + l)z w - c + z 3 ) 2 ) 2 

+ ((24 + 211 - d(z 2 + I)) 2 + ((d(z 2 + 1) + l),zi2 - c + z 4 ) 2 ) 2 

+ (23 + 213 + 1 - Z4) 2 + (24 + 214 - Mf 

+ ((*6 + 215 - *>(*3 + 2 2 )) 2 + ((6(23 + «f ) + 1)216 - a + z 5 ) 2 ) 2 

+ ((217 + 2i9 + 1 - r) 2 + (rz 20 - 25 + 2i 7 ) 2 ) 2 

+ ((zis + Z21 + 1 - r) 2 + (rz 22 - X + z 18 ) 2 ) 2 + (217 - ^is) 2 )) 2 - 



Again based on (|3.4I) . we define (with some simplifications) the polynomial £?(fc, M, a, b, r, c, d, A, X, 2, u>) as: 



(1 + A + k + 2w) 2 (2 + A + k + 2w) 2 (4w 2 + (M + 2w) 2 + 2(1 + 3w) 2 
+ ((2 + r) 2 w 2 + (1 + r + 2w) 2 ) 2 + ((A + d + w) 2 + (A + c + w + dw) 2 ) 2 
+ 2((d + 2w + dw) 2 + (c + w(2 + d + dw)) 2 ) 2 + (w 2 (2 + b + bw) 2 
+ (a + w{2 + bw(l + w))) 2 ) 2 + ((1 + r + 2w) 2 + (w + rw + X) 2 ) 2 ) 2 
+ 4(1 + 2w + 2w 2 + zf + (11 + 4A + 4k + (5 + 2A + 2k) 2 + 2z) 2 . 
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Thus, (|3.4p (and by consequence, (|3.2jl ) is equivalent to the following system of 27 Diophantine conditions 
solvable in the unknowns c, d, A, X, z, q, w, zq, . . . , z 2 4 (with k, M, a, b, r as parameters): 



P(k, M,a,b,r,c,d, A,X, z, Zq, Zi, . . . , Z24) = mod ( ^ 

\z + 1 

zo = q, 

q > z + 1 + (z + 2) z+2 {{z + 2) z+2 B{k, M, a, b, r, c, d, A, X, z, w)) w2 \ 



(3.5) ( q \ , (*i 

U + i 



Z + 1/ I JD 



We proceed as follows. Letting I be a dummy variable used only for indexing, we introduce the variables 
2/1, 2/2, ji, . . . , j 3 , /1, ■ • ■ , /24,5i, • ■ • ,524, mi, • ■ • , si, . . . , s 24 , hx,...,h 5 to represent the following equiva- 
lences: 



2/1 



<1 

z + l 



&(z + l- + ((2« + 2)« - fta(2« + 1) 2+2 - yi(2« + l) hs - h 3 ) 2 

+ ( Vl -h 4 - 2") 2 + (h 3 +h 5 + l- (2" + l)^) 2 = 0, 
P(k, M, a, b, r, c, d, A, X, z, zq, . . . Z24) = mod y\ 

O'l + h + 1 - 5i) 2 + {V\h - P( k , M, a, b, r, c, d, A, X, z, zq,... z 24 ) + ji) 2 + jf = 0, 
za = q 

(3.6) ^z o -q = 0, 

q>z + l + (z + 2) z+2 ((z + 2) z+2 B{k, M, a, b, r, c, d, A, X, z, w)) w ^ , 

<S> z + 2 + (z + 2) z+2 ((z + 2) z+2 B(k, M, a, b, r, c, d, A, X, z, w)) 10 ' 4 + y 2 - q = 0, 
For I = 1,.. . ,24, 



2/1 



(((51 + 1)2* + 1) Z! - /*((2/i + 1)2*)" - 9i) 

+ (g l+mi + l- ((y a + 1)2^) 1 ") 2 + (2/1 ^ - JO' = 0. 



From (|3.6I) . we derive our desired exponential Diophantine representation solvable in the 138 unknowns 
c, d, A,X, z,q, w, z , Z!, z 2 4,yi,V2, ji, j3, h, f24, 9i, 924,171!, ... ,m 2 4, S!, s 2 4, h 1} h 5 (with 
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k, M, a, &, r as parameters): 

((2 + 1 - h,) 2 + ((2« + 2)9 - /i 2 (2« + l) z + 2 - Wl (29 + l)' 13 - /i 3 ) 2 + (2/i - /*4 - 2 9 ) 2 + (/i 3 + fc B + 1 - (2" + l)' 11 )') 2 

+ (O'l + J2 + 1 - ?7i) 2 + (Wis - ((22 - {2A + 2k- 5) 2 - 4A - 4fc + ll) 2 + (2z - (21 + z 2 ) 2 - 3zi - z 2 ) 2 

+ {{A + k - 2 + Z23 - 21) • (A + k - 1 + Z24 - 22) • ((21 + 2 6 + 1 - 2 2 ) 2 + ((A + 27 - d) 2 + ((d + l)z 8 - c + A) 2 ) 2 

+ ((23 + z 9 - d(*i + l)) 2 + ((rf(zi + 1) + 1)210 - c + z 3 ) 2 ) 2 + ((24 + 2n - d(z 2 + l)) 2 + ((d(z 2 + 1) + 1)212 - c + 2 4 ) 2 ) 2 

+ (23 + 213 + 1 - 2 4 ) 2 + (24 + 214 - M) 2 + ((25 + 215 - b(z 3 + Z 2 )) 2 + ((6(23 + 2 2 ) + l)2l 6 - O + 2 5 ) 2 ) 2 

+ ((217 + 2i 9 + 1 - r) 2 + (rz 20 - 25 + 2i 7 ) 2 ) 2 + ((218 + 221 + 1 - r) 2 + (rz 22 -X + z 18 ) 2 ) 2 + (z 17 - z ls ) 2 )) 2 ) + ji) 2 + j 2 f 
+ (20 - q) 2 + (2 + 2 + (2 + 2)* +2 ((2 + 2) 2+2 ((l + A + k + 2w) 2 {2 + 4 + k + 2w) 2 (Aw 2 + (M + 2w) 2 + 2(1 + 3w) 2 
+ ((2 + rfw 2 + (1 + r + 2w) 2 ) 2 + {{A + d + wf + (A + c + w + dw) 2 ) 2 + 2((d + 2w + dw) 2 + (c + w(2 + d + dw)) 2 ) 2 
+ (w 2 {2 + b + bw) 2 + (a + w(2 + bw(l + w))) 2 ) 2 + ((1 + r + 2w) 2 + (w + rw + X) 2 ) 2 ) 2 + 4(1 + 2w + 2w 2 + z) 2 

24 

+ (11 + AA + 4fc + (5 + 2A + 2k) 2 + 2z) 2 )) w2i + y 2 - qf + £ ( (((y a + 1)2* + l) 2 ' - / ; ((yi + 1)2 Z T ™ 9if 

1=1 

+ (g l+ m l + l- ((yi + 1)2")") 2 + (yiSi - /0 2 ) 2 = 0. 

The use of summation notation in the above representation is permissible, and in fact is more informative 
than expanding it. For instance, Keijo Ruohonen [16] . uses summation notation in his 79 variable Diophan- 
tine representation of Fermat's Last Theorem in order to keep it under a page. In fact, Davis, Matiyasevich, 
and Robinson [3j pg. 332] acknowledge this as a writeable representation. But it is not a problem either if 
we were to expand the sum, as the entire representation would cover exactly half a page. 

Now, in order to obtain a Diophantine representation of PH 2 , it would not be sufficient to simply replace 
every exponential function with its Diophantine representation (say, for the sake of example, even a rep- 
resentation that is economical with respect to the number of variables, namely, five, as presented in [11] V 
since this would result in a representation of 638 variables. Instead, if one uses the 10 variable Diophantine 
representation of the binomial coefficient obtained via Theorem [7] (and replace all the exponential functions 
of the sufficiently strong inequality for q in ()3.6|) with the five variable representation) , then we will obtain a 
Diophantine representation of PH 2 in 347 variables (though if one is concerned with writeability and not just 
the theoretical minimization of variables, then the 18 variable version should be used instead). We should 
also mention that for the specific case of PH 2 , the exponential Diophantine and Diophantine equations ob- 
tained could further be used via some combinatorial tricks in Bovykin and De Smet's original intermediate 
representation of PH 2 presented in (|3.1j) . which Bovykin and De Smet are apt to point out. However, despite 
the naivete of the starting representation in (|3.1[) prior to the elimination of the bounded universal quantifier, 
the fact that we have still obtained compact exponential Diophantine and Diophantine representations of 
PH 2 points to importance of the conservation of variables due to the techniques of §2. To compare, had we 
directly used the method of Davis, Putnam, and Robinson [2] (namely, (|2.2[) and (|2.3jl ) without any of the 
techniques in §2 to eliminate the bounded universal quantifier in (|3.4j) . then a rough estimate would yield an 
exponential Diophantine representation of 233 variables and, under the assumption of a five variable repre- 
sentation of the exponential function, a 1055 variable Diophantine representation. Thus, we have managed 
to conserve 95 variables and 708 variables for each case, respectively. 

3.2. Proof of Theorem [4l Goodstein's Theorem [4] is the following: 

Theorem 9 (Goodstein's Theorem). Given any non- decreasing function p r , po>2, a number no, and the 
function n r defined as: 

n r+1 = S$r +1 (n r ) - 1, 

then 3r. n r = 0, where (n r ) is the operation of putting n r in hereditary base-p r notation and then 
replacing every occurence of p r in this representation with p r +i- 
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We will instead use the following simpler notation (due to Kir by and Paris 7J): Let G n (m) be the number 
produced by replacing every n in the hereditary base-n representation of m by n + 1 and subtracting 1. So, 
the Goodstein sequence for to starting at 2 is: 

too = to, to-i == G 2 (m ), to 2 = G 3 (toi), to 3 = G 4 (to 2 ), . . . 

Hence, rrii can be viewed as the cipher of the base i + a positional code (to^, i + a, I + 1) of the sequence 
(co, . . . ,c;}. Thus, the statement of Goodstein's theorem can be rewritten as (where m > 1 and a > 1 are 
taken as parameters): 

3k Vi < k (rrii + i — Gi+ a (mi) A too = to A TOfc = 0). 

In other words, the Goodstein sequence for m starting at a eventually terminates. Hence, the bulk of the 
work will be focused on representing the action of Gi+ a (TOi) as a prefixed polynomial expression. 
For any rrii , the base i + a representation of rrii is 

to 2 = c t (i + a) 1 + q_i(i + a) 1 ^ 1 + . . . + ci(i + a) + c . 

Then Gi+ a (jn-i) can be defined in terms of a recursive function / (which puts rrii in its hereditary base i + a 
representation and replaces every occurence of i + a with i + 1 + a) as such: 

G i+a W = r'' + "(« + l + «)-l, 

where 

(3.7) f mi ' i+a (i+l + a) =Y J Cj(* + ^ + a) f3 ' l+a{l+1+a \ 

3=0 

where for each Cj (0 < j < I) has the following Diophantine representation, namely, 

Cj = Elem(TOj, i + a,j). 

Note that by definition, 

f°- l+a {i + l + a) :=0. 

So the next step is to construct a prefixed polynomial expression for /. The first step is to view /, the 
process of putting a number in hereditary base notation, in terms of levels that each terminate at the highest 
power of the level. Note that the highest power I of the base i + a representation of a natural number n can 
be given as the following Diophantine representation: 

I = HP(n) ((i + a) l+1 > n A n > (i + a) 1 A n ^ 0) V (n = A I = 0) . 

I suspect though that there must be a number-theoretic function that already does this, so I do not need 
as many variables in the representation of I = HP(n) (however, for now, this will have to suffice). With 
this definition in hand, we can proceed to define ReplaceBase(rt), which returns the natural number which 
results from taking the base i + a representation of n and replacing every occurence of i + a with i + 1 + a: 

n' = ReplaceBase(n) <^> 

31 Vk<l 3c [I = HP(n) A a = Elem(n, i + a, 0) A a k+1 = a k + c ■ (i + a + l) k+1 

A c = Elcm(?i, i + a, k + I ) A a; = n ]. 

Note that the sequence ao, . . . , a/ will be Godel coded (by the pair b, d) so the above definition of ReplaceBase 
formally becomes 

n = ReplaceBase(n) <^> 
3lbd \fk<l3c[l = HP(n) A rem(6, 1 + d) = Elem(n, i + a, 0) 
A rem(6, 1 + (k + 2)d) = rem(6, 1 + (k + l)d) + c • (i + a+ l) k+1 
A c = Elcm(n, i + a, k + 1 ) A rcm(6, 1 + (I + l)d) = n']. 
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For clarity of exposition, I will leave sequences be and not replace them with the Godel code of their elements 
via the remainder function. So, / can be viewed as the following sequence, So, Si, . . . , sl where 

s = RB(m i ,(0,...,HP(m i ))) 

Sl = RB(m l; (RB(0, (0, . . . , HP(0)», . . . , RB(HPK), (0, . . . , HP(HP(m,-))»» 

(3.8) 

s L = m+i + 1. 

Namely, sq is the ReplaceBase operation applied to m% and (0, . . . , HP(m,)) are the exponents of sq. Next, si 
is the natural number that results from taking sq and applying the ReplaceBase operation to its exponents. 
S2 is the natural number that results from taking sq and applying the ReplaceBase operation to its exponents. 
And so forth, until we terminate at Sl, which gives us the natural number m i+1 + 1 (since when we subtract 
1 we get rrii+i). First, we can determine L, more exactly as: 

L = level(mj) 

Vfc < L[ao — rtii A a k+ i = HP(afc) A < i + a A cll-i > i + a]. 

It is also important to be able to access each of the exponents of any s n in the sequence defined in (|3.8[) 
because then we can apply the ReplaceBase operation on them. Thus, the property ll p is the fc-th exponent 
of s„" is Diophantine: 

P = Ex Pfc( s «) 
3cy [c(i + 1 + a) p + y — s n A c = Elem(mi, i + a, k)]. 

With these two representations, the function f ni ' 1+a (i + l + a) (as defined in (|3.7[l ) can be represented as the 
following prefixed polynomial (note that we use the following three sequences and omit the explicit use of the 
remainder function to represent their elements for clarity of presentation: z , . . . , Zl, Sq, ■ ■ ■ , sl, t , . . . , tj): 

3L Vn < L 31 Vfc < / [zq — mi A z n +i = HP(z„) A zl < i + a A zl-i > i + a 

A sq = ReplaceBase(mi) A I = z\ A To — c • (i + 1 + a) d Ac = Elem(mi, i + a, 0) 

A d = ReplaceBase(d') Ad' = Exp (s„) A r k+1 = r fe + / • (i + 1 + a) d 

A f = Elem(mj, i + a, k + 1) A d'" = ReplaceBase(ci") A d" = Exp fc+1 (s„) A n — s n+ i A s L = m t+ i + I]. 

Thus, the statement of Goodstein's theorem is: 

3r \/i < r 3L Vn < L 31 Vfc < I [zq — rrii A z n+ \ = HP(z„) A zl < i + a A zl-i > i + a 

A s a = ReplaceBase(rrii) A t = c • (i + 1 + a) d Ac — Elem(mi,i + a, 0) 

Ad = ReplaceBase(d') Ad' = Exp (s„) A r k+1 = r k + f ■ (i + 1 + a) d "' 

A f = Elem(mi, i + a, k + 1) A d'" = ReplaceBase(d") A d" — Exp fc+1 (s„) A n — s n+ i 

A sl — TOi+i + 1 A mo = m A m r = 0]. 

Expanding this, we get the following 181 variable exponential prefixed polynomial representation: 
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3rbidi Vt < r 3viVi 4 vi 5 viQhih2h 3 h 4 w 3 w 4 piP2L Vn < L 3v2V 3 v 4 v^VqVsVi2Vi 3 WiW2W^WqW 7 Wsw 9 WiqWiiWi2Wi 3 Wi 4 Wi^Wiq 
w\ 7 jihhhhh e i e 2Uiu 2 u 3 Vu 5 < tii 3u 4 uq . . .u 30 l Vfc < I 3v 7 cdd'v 9 v w fd"d'"h^hQh 7 hsh 9 hiohiihi 2 hi 3 hi 4 V/115 < /112 

3ftl6 ■ • • 7141.77 ■ • • hi Vj25 < .722 3j 2 6 • • • J5l9l • • • 99 

[((«i + u> 3 - di(i + l)) 2 + ((di(i + 1) + l)w 4 - 61 + wi) 2 ) 2 + ((«2 + u> 5 - d 3 (ra + 2)) 2 + ((d 3 (n + 2) + l)io 6 - 6 3 + « 2 ) 2 ) 2 

+ ((«s + «>io - d 3 (n + I)) 2 + ((ds(n + 1) + ljttfn - 63 + «3) 2 ) 2 + ((«i +101- d 3 ) 2 + ((d 3 + l)u>2 - & 3 + «i) 2 ) 2 

+ ((((«3 +w 7 + l-(i + ar +1 f + ((t + aD 2 + w s - v 3 ) 2 + (4 -w 9 - l) 2 ) • (« 2 + vl)) 2 + ((v 4 + u, 12 - d 3 (L + l)) 2 

+ ((d 3 (L + 1) + l)w 13 - 63 + w 4 ) 2 ) 2 + (v 4 + ttfi4 + 1 - » - a) 2 + ((«6 + wis - d 3 L) 2 + ((d 3 L + l)to 16 - b 3 + v 5 ) 2 ) 2 

+ (i + a + wn + 1 - V5) 2 + ((ve + ji - d 2 ) 2 + ((d 2 + l)j 2 - b 2 + v e ) 2 ) 2 + ((u e + u w - u 3 ) 2 + ((u 3 + l)uu - u 2 + u e ) 2 ) 2 

+ ((uw + l) 2 + (U12 + u 6 (i + a) Uli + U13 - wi) 2 + K + wig + 1 - (* + a)) 2 + («is + «ie + 1 - (» + a)" 14 ) 2 ) 2 

+ ((u 7 + U17 - (U 8 + 2)u 3 ) 2 + (((u 8 + 2)« 3 + l)ui 8 - ti 2 + M 7 ) 2 ) 2 + (K + «19 - («6 + 1)" 3 ) 2 + (((u B + 1)« 3 + 1)«20 

-u 2 + u s ) 2 ) 2 + (tig + u 4 (i + a + 1)" 5+1 - w 7 ) 2 + ((u 5 - U23) 2 + (W2i(i + a)" 5+1 + « 4 (« + a)" 23 + "22 - W4) 2 

+ (m 4 + «24 + 1 - (i + a)) 2 + (?i 22 + M25 + 1 - (i + a)" 23 ) 2 ) 2 + ((tig + U 26 - (til + l)ti 3 ) 2 + ((("I + l)t*3 + 1)^27 

- U2 + M 9 ) 2 ) 2 + (tig - M 6 ) 2 + ((((«! + M 28 + 1 - (t + «)" 1 + 1 ) 2 + ((» + «) M1 ) 2 + «29 - «l) 2 + («? - «30 - I) 2 ) ' («? + " 2 )) 2 ) 2 

+ (K + ^5 - d 4 ) 2 + ((d 4 + l)/i 6 - "2 + «?) 2 ) 2 + («r - c • (i + 1 + aff + ((/i 9 + I) 2 + (h 7 + c(i + a) h ° + h s - Vl ) 2 

+ {c + h w + l-(i + a)) 2 + (hs + /in + 1 - (* + a) hs ) 2 ) 2 + ((h 17 + h 21 - h 14 ) 2 + {{h 14 + l)h 22 - h 13 + h 17 ) 2 ) 2 

+ ((/i26 + I) 2 + (h23 + h 17 (i + af™ + h 2 ± - d') 2 + (h 17 + h 26 + l-(i + a)) 2 + (h 24 + h 27 + 1 - (i + a)' 125 ) 2 ) 2 

+ ((his + h 2S - h 14 (h 15 + 2)) 2 + ((h 14 (h 15 + 2) + l)h 29 - h 13 + h ls ) 2 ) 2 + {(h w + h 3 o — hi4(hi5 + 1)) 

+ ((hu(h 15 + 1) + l)h 31 - h 13 + h 19 ) 2 ) 2 + (h 19 + h 16 (i + a + - h w ) 2 + ((/lis - /i 34 ) 2 + (/i32(* + a) fel5+1 

+ /i 16 (i + a)' 134 + /i 33 - d') 2 + {hits + h 35 + 1 - (* + a)) 2 + (h 33 + h 36 + 1 - (i + a)* 134 ) 2 ) 2 + ((h 20 + h 37 - h u {h 12 + l)) 2 

+ ((h 14 (h 12 + 1) + l)/i38 - /lis + ^2o) 2 ) 2 + (^20 - d) 2 + ((((d' + /i 39 + 1 - (« + a)' ll2 + 1 ) 2 + ((i + a)' 112 ) 2 + /i 40 - d') 2 

+ {{d'f - h 41 - l) 2 ) • ((d') 2 + h\ 2 )) 2 ) 2 + ((vs + ei - (n + l)d 2 ) 2 + ((d 2 (n + 1) + l)e 2 - fe 2 + m 8 ) 2 ) 2 + ((j'r(i + 1 + af 

+ js - vs) 2 + ((.7ii + I) 2 + (.79 +jr(i + a)' 11 + j w - Vl ) 2 + (j 7 + j 12 + 1 - (i + a)) 2 + (j w + j 13 + 1 - (i + a)^) 2 ) 2 ) 2 

+ ((«9 + - (k + l)d 4 ) 2 + ((d 4 (fc + 1) + l)ji 5 - fe 4 + w 9 ) 2 ) 2 + ((wio + jie - (fc + 2)d 4 ) 2 + ((d 4 (k + 2) + l)ji 7 - 6 4 

+ «io) 2 ) 2 + («9 +/•(* + 1 + af" - v w ) 2 + ((jig - fc) 2 + (ji 7 (i + a) k+1 + f(i + a)^ + j 18 - «i) 2 + (f + jm + 1 - (i + «)) 2 

+ (.718 + .721 + 1 - (* + a) Jl9 ) 2 ) 2 + (C?27 + 331 - hif + ((.724 + 1).732 - J23 + j 27 ) 2 ) 2 + ((j 3 5 + I) 2 + (.733 + ,72 7 (« + a)^ + j 34 

- d") 2 + (.727 + .736 + 1 - (* + a))' + (.734 + J37 + 1 - (t + a) j35 ) 2 ) 2 + ((.728 + J38 - .724^25 + 2)) 2 + ((.724 (.725 + 2) + l)j 3 9 

- .723 + ,72 8 ) 2 ) 2 + ((.729 + .740 - hi(.]25 + l)f + {{hiijV, + 1) + l)j41 - J23 + J2 9 ) 2 ) 2 + (.729 + .726 + a + l)^ 1 - ,72 8 ) 2 

+ (C726 - .?44) 2 + (Mi + a)' ?25+1 + hS + a) Hi + J 43 - d") 2 + (j 26 + .745 + 1 - (* + a)) 2 + (j 43 + j 46 + 1 - (i + a)^ ) 2 ) 2 

+ ((h0 + .747 - .724(722 + I)) 2 + ((.724(722 + 1) + 1).748 - J23 + .73o) 2 ) 2 + (ho - d'"f + ((((d" + j\ 9 + 1 - (i + a)-' 22 + 1 ) 2 

+ ((t + ap 2 ) 2 + j 50 - d") 2 + ((d") 2 - hi - I) 2 ) • ((d") 2 + .7 2 2 2 )) 2 ) 2 + ((«12 +qs~(l + l)d 4 ) 2 + ((d 4 (l + 1) + 1)59 
-b 4 + v 12 ) 2 ) 2 + ((qi(i + 1 + a) d " + q 2 - v s ) 2 + ((q 5 - k) 2 + (q 3 (i + a) k+1 + qi(i + a) q5 + q 4 - Vl ) 2 + (qi + q 6 + 1 

- (* + a)) 2 + (q 4 + q 7 + l-(i + a)^) 2 ) 2 ) 2 + ((v 12 + j 3 - d 2 (n + 2)) 2 + ((d 2 (n + 2) + l)j 4 - b 2 + v 12 ) 2 ) 2 + ((u 13 +j 5 
-d 2 (L + l)) 2 + ((d 2 (L + 1) + l)j 6 -b 2 + v 13 ) 2 ) 2 + (vi 3 - v u - l) 2 + ((v m +Pl -di(i + 2)) 2 + ((di(i + 2) + l)p 2 
-h+ v u ) 2 ) 2 + (v 14 + 1 - w 13 ) 2 + ((u 1B + hi - di) 2 + ((di + l)h 2 - h + vi 5 ) 2 ) 2 + (v 15 - m) 2 + ((v 16 + h 3 - di(r + l)) 2 
+ ((di(r + 1) + l)/i 4 -h + v 16 ) 2 + v 2 6 ) 2 = 0]. 
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